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Abstract. High-resolution Ha monitoring of Cyg X-1, HD226868 was carried out during 1996-2002 and the 
resultant spectra analysed in conjunction with 1.5-12 keV X-ray monitoring. We demonstrate that the Ha line- 
profiles have complex variability on different timescales, controlled in particular by the orbital period and the 
focused wind model of mass loss. We find that long-term variability of the mass loss by the supergiant and short- 
term variability due to clumpy structure of the stellar wind dominate during the low/hard X-ray state and that 
X-ray photoionization has a relatively small influence on the line-profile shape and EW variability. During the 
high/soft X-ray state and flaring the effect of photoionization the line-proflle and EW of Ha increases but is 
still unable to describe the loose anti-correlation between EW and the low energy X-ray emission. We propose 
that variability of the mass loss by the supergiant can change wind velocities in the Stromgren zone around the 
accretion disc of the secondary, leading to an increase in accretion rate through the disc. 
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1. Introduction 

Cyg X-1 has been one of the most well-studied X-ray 
sources since its discovery in 1965 (Bowyer et al. 1965) 
and identification with the 09.7Iab supergiant HD226868 
(Bolton 1972; Webster & Murdin 1972). It provided the 
first observational evidence for the existence of black holes 
when it was found to be a spectroscopic binary with an 
orbital period of 5.6 days. The ephemeris was revised a 
number of times and a recent estimation based on spec- 
troscopic and photometric orbital variability is presented 
by Brocksopp et al. (1999a). The masses of the two com- 
ponents were calculated as 17.5 Mq for the supergiant and 
10.1 M0 for the black hole, with adopted system inclina- 
tion i=35°(Herrero et al. 1995). 

Cyg X-1 is a detached double system but the optical 
component is very close to filling its Roche lobe. Typically 
of high mass stars, HD22686 is thought to emit a stellar 
wind, which appears to be attracted towards the black 
hole as in the focused wind model of Friend & Castor 
(1982; Gies & Bolton 1986; Sowers et al. 1998). The opti- 
cal spectrum of the system has only two well-pronounced 
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emission lines: Ha and HeiiA4686, but high-resolution 
spectroscopy detects line-asymmetry in most of the hy- 
drogen and red Hei lines (Canalizo et al. 1995) which 
thus cannot be used for precise determination of the or- 
bital elements. Gies & Bolton (1986) and Ninkov et al. 
(1987) found that the Heii A4686 emission component has 
a radial velocity curve that is best interpreted as emission 
originating between the stars. The Ha profile appears to 
have at least two components; P Cygni emission moving 
with the orbital motion of the primary star and a sec- 
ond emisison component which follows the quite different 
radial velocity curve of the He ii A4686 emission line and 
which is formed in the focused wind flow between the stars 
(Ninkov et al. 1987; Sowers et al. 1998). 

Cyg X-1 demonstrates sporadic flaring at soft X-ray 
energies (2-12 keV) on different timescales, from days to 
months and years (Brocksopp et al. 1999b). It is believed 
that the source of that variability is located near the black 
hole companion and is connected with the accretion disc. 
A supergiant wind is also variable on different timescales 
from days to years (Kaufer et al. 1996; Rivinius et al. 1997) 
and the variability of the emission lines and the soft X-rays 
should be somehow connected because they both depend 
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Table 1. Equivalent widths and orbital phase values of the Ha line-profiles. Spectra marked by an asterisk have been 
plotted in Figs. 1 and 2. 
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on mass transfer via the wind. Discovery of any correla- 
tion between the stellar wind and the X-ray variability 
is important for understanding the physical processes of 
mass transfer through the disc and onto the black hole. 

Long-term monitoring of variability of the stellar wind 
from the supcrgiant of Cyg X-1 has not bc;en previously 
attempted. Nonetheless some orbit-independent variabil- 
ity of the equivalent widths of the HeiiA4686 emission 
line was found by Ninkov et al. (1987) and of the emission 
component of Ha by Voloshina et al. (1997). The latter 
also discovered a decrease of EW accompanying the period 
of high/soft state behaviour in 1996 (see also Brocksopp 
et al. 1999b). Thus we have monitored the Ha line at 
the Crimean Astrophysical Observatory (CrAO) since the 
1996 X-ray spectral state change, when a strong increase 
in soft X-rays was found by the Rossi X-ray Timing 
Explorer (RXTE); our observations were supported by 



long-term photometric monitoring provided by Sternberg 
Astronomical Institute of Moscow. 



2. Spectroscopic observations and data reduction 

The spectroscopic monitoring of Cyg X-1 took place over 
seven years from 1996 to 2002, during which time we ob- 
tained a total of 142 spectra. We used the Coude focus 
of the 2.6-m telescope of the CrAO and the detector was 
an Electronix CCD array (1024x260 pixels). All observa- 
tions were performed in the second order of a diffraction 
grating with a reciprocal dispersion 3 A/mm. The spectral 
resolution was about 35000 and during each observation 
we obtained a 60 A spectrum. The exposure time ranged 
from 30 to 45 minutes, depending on weather conditions 
resulting in a signal-to-noise ratio of 100 or better. In most 
cases two or three spectra were obtained during each night 
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Fig. 1. Selected Ha profiles for orbital phases 0.00-0.50, determined according to the cphcmcris of Brocksopp ct al. 
(1999). Each box contains three sample spectra, chosen at random in terms of continuum or X-ray activity but so as 
to include extremely strong, weak and intermediate intensity examples of the Ha emission line in each case (plotted 
spectra are indicated in Table 1). Spectra drawn with a thick line are those obtained during periods of X-ray activity. 
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Fig. 2. Selected Ha profiles for orbital phases 0.50-1.00, chosen as in Fig. 1 
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and then added for improved signal-to-noise and cosmic 
ray subtraction. 

The spectra were reduced using standard techniques 
for dark-current subtraction and flat-field division. The 
subsequent reduction included removal of spurious spikes 
produced by cosmic-rays, sky-background subtraction, 
water vapor line subtraction using non-emission fast- 
rotating early type stars and filtering of the spectra by 
the moving- average method using three points. For the 
wavelength calibration we used ThAr comparison spec- 
tra obtained immediately before and after the target ob- 
servation. The accuracy of the calibration was 1.5 km/s 
or better. All spectra were corrected to the baricentre of 
the solar system, both in wavelength and time of mid- 
exposure. The continuum normalization was provided by 
fitting of a fifth degree polynomial. 

As a final step, equivalent widths (EW) of each spec- 
trum were measured. Estimated errors of EW measure- 
ments were better than 0.1 A. The results are presented in 
Table 1. 

During the period of our observations the system was 
monitored by RXTE. In order to compare variability of 
Ha with the X-ray variability, we use the All Sky Monitor 
(ASM) pubhc archive data from the web. The ASM ob- 
servations were made in the three energy bands: A (1.5 3 
keV), B (3-5 keV) and C (5-12 keV). A detailed descrip- 
tion of the ASM, including calibration and reduction, is 
published in Levine et al. (1996). The daily mean data 
were converted from counts to energy units (keVcm~^s~^) 
using transformation formulae obtained from Zdziarski et 
al. (2002). 

3. Ha line-profile variability 

Previous study of the emission spectrum of Cyg X-1 shows 
that the Ha line-profile has strong variability on the 5.6- 
day orbital period (Gics & Bolton, 1986; Ninkov et al. 
1987; Sowers et al. 1998). However our observations show 
high variability of the emission line intensity that is phase- 
independent and which was not noticed in the previous 
observations (which took place over a shorter time range). 
Since mass loss by the optical component is the only source 
of accretion, any instability of the stellar wind will lead 
to a variable mass accretion though the response of the 
accretion disc and that will have some impact on the X- 
ray emission. 

Such complex variability makes analysis of the line- 
profile variations much more difficult to interpret because 
the shape and strength of the Ha emission profile depends 
on a number of factors: 

i) phase of the orbital period (wind is noticeably coUi- 
mated in the direction of the secondary) 

ii) variability of the stellar wind density due to variable 
mass loss rate by the supergiant 

iii) clumpy structure of the stellar wind providing discrete 
absorption components (DAC) 



iv) time-variable and phase-locked photoionization of the 
stellar wind by the soft X-ray radiation from the ac- 
cretion disc of the secondary component 

In order to overcome this we have divided our spectra 
into orbital phase groups spanning 0.05 in phase and anal- 
ysed each group separately according to processes that 
may have influenced them. Figs. 1 and 2 show three sam- 
ple line-profiles for each group, chosen at random in terms 
of continuum or X-ray activity but so as to include ex- 
tremely strong, weak and intermediate intensity examples 
of the Ha emission line in each case (plotted spectra are 
indicated in Table 1). Spectra drawn with a thick line are 
those obtained during periods of X-ray activity (i.e. inten- 
sity of the 1.5-3 keV X-rays in the exceed 4 keVcm~^s~^) 
and it appears that the X-ray variabilty has little effect 
on the profiles (see Sect. 4). 

Figs. 1 and 2 show that the orbital variability of the 
line profile is mostly in agreement with the focused wind 
model; there are two emission components as demon- 
strated by Sowers et al. (1998). One forms a relatively 
stable PCyg structure, which is produced by the stellar 
wind from the hot supergiant and moves around the spec- 
trum in the phase of the orbital period. The other emission 
component originates between the stars in a focused wind 
flow from the supergiant to the unseen companion. The 
orbital movement of this component of the line-profile is 
different since it appears in the coUimated stellar wind 
close to the LI point, such that it is blue-shifted dur- 
ing the phase range 0.10-0.60 and merged with the P Cyg 
emission component at other times. 

Irregular variability of the line-profiles is clearly seen 
during all phase intervals; sometimes (phases: 0.25-0.35, 
0.50-0.55, 0.70-0.75) there is an increase or decrease in 
the intensity of the whole line but for the majority of 
phase intervals irregular variation of the blue part of the 
profile seems to dominate. These variations can be con- 
nected with the variability of either the collimated wind 
component or with the absorption part of the P Cyg line- 
profile from the stellar wind of the supergiant. However 
since both emission components reflect variable mass loss 
from the optical star, it is most probable that they are 
connected with one another. 

Fig. 3 shows long-term light curves for the RXTE in- 
tensity (Band A), hardness ratio (HRl = Bands B/A) 
and Ha EW. As with the line profiles it is clear that 
there is complex variability in the EW measurements, on 
both short (~days e.g. MJD 51392-51398) and long (~few 
weeks e.g. MJD 51980-52200) timescales but the sampling 
is often insufficient to be conclusive. However, it is certain 
that this variability is not connected with the orbital pe- 
riod and instead seems to reflect some other mechanism(s) 
within the system. 

4. Connection with the soft X-ray variability 

The ASM data plotted in Fig. 3 demonstrate that after the 
high/soft state period (i.e. the major outburst in the softer 
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Fig. 3. Long-term variability of: a)X-ray ASM intensity in the 1.5-3 keV band (keVcm ^); b) hardness ratio (HRl) 
of the soft X-ray radiation, (3-5 kcV)/(1.5-3 kcV); c) EW of Ha (A). Filled circles on Figs. 3a & 3b indicated X-ray 
measurements obtained simultaneous with Ha observations. 



X-ray bands) in 1996 (MJD 50100-50400) the subsequent 
flares increased in both intensity and softness from year to 
year, reaching a maximum in 2002 with another high/soft 
state period. These flares have a semi-regular nature with 
characteristic occurrence time of ~400 days. 



The variability of the EW and ASM X-ray flux are of- 
ten loosely correlated. For example the most active X-ray 
periods MJD 50200-50350 and 52100-52500 were both ac- 
companied by a decrease in EW; similarly the more short- 
lived flares during MJD 51850-51900 and 51440-51480. 



A.E.Tarasov et al.: Variability of the Ha emission of Cygnus X-1 



7 



15 



10 



..« • 



I ! I L 



m 
m 



CD 



a 

CO 

m 




1 1 

• 


1 

c 


• 

• • 


• — 






1 1 


1 



EW 

Fig. 4. Soft X-ray intensities plotted against the EW of 
Ha for each of the three X-ray energy bands. X-ray inten- 
sities. Characteristic behaviour can be seen for some (but 
not all) ranges of EW for Band A and not at all for Band 
C. 



However a decrease in EW is not proportional to the am- 
plitude of X-ray flares and is not always associated with 
episodes of X-ray activity. Furthermore the EW of Ha 
decreased monotonically from —1.1 A on MJD 51397 to 
0.3 A on MJD 51456 (Fig. 3b) prior to the X-ray flare, 
reaching a minimum simultaneously with the small X-ray 
peak. 

To understand the situation better, we combined the 
EW of Ha with the three different ASM spectral bands. 
Fig. 4 shows that the relationship between the EW and 
the softer X-ray emission depends strongly on the nature 
of the spectral state; X-ray flares appear only if the EW 
drops below —0.4 A. The softest X-ray spectral band (A) 
is most sensitive to the Ha variability such that a de- 
crease in EW will probably, although not necessarily, be 
associated with an increase in X-ray flux as seen in the 
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Fig. 5. Relationship between the hardness of the low- 
energy X-ray spectrum and EW of Ha; a) HRl - (3-5 
keV)/(1.5-3 keV) and b) HR2 - (5-12 keV)/(3-5 keV). 
This apparent correlation is clearly more signiflcant for 
HRl (top) than HR2 (bottom). 



lightcurves; on the contrary the hardest X-ray band (C) 

appears insensitive to the EW. 

The hardness of the X-ray spectrum can be a useful 
parameter for understanding the relation between X-ray 
emission and emission in the Ha line - the EW and X- 
ray hardness ratio appear to be correlated for both X- 
ray colours (Fig. 5), although that of HRl (i.e. the softer 
colour) is more significant. The scattering of the data in 
Fig. 5 is significantly larger than the estimated errors of 
the observations, both in EW and X-ray and thus some 
additional factors must play an important role in the con- 
nection between optical and X-ray variability. 

Comparing Figs. 4 and 5 it is possible to conclude that 
the loose correlation between the Ha and X-ray emission 
is more significant for the X-ray spectral hardness than 
for the intensity. 

5. Discussion 

A number of sources could be responsible for the pres- 
ence of some relationship between the variability of the 
Ha EW and low-energy X-ray radiation; these include (i) 
photoionization of the neutral gas by the soft X-ray radi- 
ation, (ii) X-ray flares and state transitions of the accre- 
tion disc due to variable mass transfer via the stellar wind 
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Fig. 6. All our Ha line profiles obtained during phase intervals 0.95-1.10 (a) and 0.45-0.55 (b). Profiles marked by an 
asterisk were obtained during the time of X-ray activity, when intensity in the band A (1.5-3.0 keV) were more than 



and/or (iii) variability of X-ray radiation due to variable 
column density of neutral gas in the line of sight. 

5.1. Photoionization vs. mass-loss variability 

Unlike single hot supergiant stars, Cyg X-1 has an addi- 
tional source of ionisation in the form of X-ray radiation 
from the accretion disc of the secondary. However, due to 
orbital motion, it is phase locked with the orbital period 
and can be identified relatively easily. This is illustrated 
in Fig. 6 where we display all spectra in the orbital phase 
intervals 0.95-1.10 (left) and 0.45-0.55 (right). During the 
first of these phase intervals the black hole is situated be- 
hind the optical component, with the red component of 
the Ha profile formed in the part of the envelope which 
is directed towards the X-ray source (resulting in a su- 
perposition of the envelope and focussed wind emission 
components) . During the second phase interval the star is 
behind the black hole and the red emission is no longer 
under the influence of the X-ray source. 

Those spectra obtained during periods of X-ray activ- 
ity have been indicated in Fig. 6 by an asterisk and show 
clearly that the effects of photoionisation by the X-ray 
emission do not dominate the shape of the emission line. 
However, comparison of the red wing of the red emission 
peak at different orbital phases shows that photoionisa- 
tion does have some influence on the line-profile. At phase 
0.0 the red wing has negligible extension and the photo- 
spheric line profile is often seen; this can be attributed to 
photoionization by the X-rays since the red wing of the 



emission forms bewteen the star and the black hole. At 
phase 0.5 the red wing forms on the far side of the star 
where it is unaffected by the X-ray radiation and can be 
extended up to 700 km/s. During both phase intervals 
plotted in Fig. 6 the aperiodic part of the wind variability 
has a strong influence over the shape of the Ha line. Often 
this sporadic inhomogeneity of the wind or long-term vari- 
ability of mass loss is so high that it prevails over the 
(focused wind model controlled) orbital line-profile vari- 
ability. This can occur at any phase interval of the orbital 
period and independently of X-ray activity. 

In the top panel of Fig. 7 we plot the EW values against 
orbital phase, using filled circles for the data obtained 
during periods of X-ray activity (intensities in the 1.5- 
3 keV band exceed 4 keVcm~^s~^). While the quiescent 
(and combined) data show only scatter, the scatter for the 
filled circles (i.e. epochs of X-ray flaring) is much more 
constrained. 

The presence of the orbital variability is best seen 
in the EW of the red emission peak alone. The highly 
scattered wave in Fig. 7b (non-flaring episodes) is in full 
agreement with the focused wind model and increasing of 
EW at phase 0.6-0.9 can be seen as the collimated wind 
component merges with the P Cyg wind component (see 
Fig. 2). The photoionization shown in Fig. 6 is relatively 
small in amplitude and the large scattering of the points in 
Fig. 7b can be best-described by the short- and long-term 
variability of the stellar wind. Fig. 7c (flaring episodes) 
shows that photoionisation dominates during phases 0.6- 
1.3 when the X-ray source is located close to where the 
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Fig. 7. (a) Variability of the equivalent width of the Ha 
line with the phase of orbital period. Open circles - EWs 
obtained outside X-ray activity; filled circles - EWs ob- 
tained during X-ray activity, (intensities of the soft X-ray 
in the 1.5-3 keV band exceed 4 keVcm~^s~^). Typical 
error of EW measurements are shown for one of points, 
(b) Similar plot showing the EW of the red emission peak 
only, during non-flaring episodes, (c) Similar plot show- 
ing the EW of the red emission peak only, during flaring 
episodes. 



red emisison peak is forming. Both curves on Figs. 7b & 
7c are statistycally significant the full amplitudes of EW 
variability are exceed more than ten times standard devi- 
ations of the data. This is in agreement with the orbital 
variability of the resonance UV lines studied by Thevcs ct 
al. (1980) and van Loon et al. (2001). The resonance UV 
lines form in a more extended region (which includes the 
Roche lobe of the secondary) than Ha and are significantly 
influenced by photoionization even during quiescence. 

The complex relationship between the X-rays and Ha 
EW that was seen in Figs. 4 and 5 cannot be explained 
purely in terms of time- variable photoionisation. Plotting 
EWs of the red emission peak against the 1.5-3 keV X- 
rays for phase intervals 0.0±0.2 (filled circles) and 0.5±0.2 
(open circles) shows a similar plot to Fig. 4 for each phase 
interval (Fig. 8). If photoionization dominated then we 
would expect to see correlated behaviour in each phase 
interval; instead a large scatter and independent variations 
are still present. 

The decay of EW prior to the X-ray flare of MJD 51456 
would suggest that (at least in some instances) a decrease 
in mass loss can trigger an X-ray flare. A reduction in 
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Fig. 8. Soft X-ray intensities (1.5-3 keV energy band) 

plotted against the EW of Ha red emission peak. Filled 
circles - phase interval of the orbital period 0.0 ±0.2; open 
circles phase interval 0.5 ± 0.2. 



mass-loss would increase the size of the Stromgren (i.e. 
photoionisation) zone surrounding the X-ray source (van 
Loon et al. 2001). This would increase the accretion rate 
through the disc and result in X-ray activity, thus increas- 
ing the size of the Stromgren zone further. In the case of 
the 2001-2002 high/soft state of Cyg X-1, once the EW 

return to —1.5 2.0 A the Stromgren zone would be 

compressed to its normal size, as was the case after the 
1996 high/soft state . However, this would appear not to 
be the case for all flares and optical observations with 
improved time resolution are needed to investigate this 
further. 

5.2. Variability of the column density toward secondary 

It is well known that the winds of early type stars are 
driven by resonance line scattering of radiation in the stel- 
lar envelope. This radiative driving mechanism is inher- 
ently unstable. Small perturbations of an otherwise mono- 
tonic wind velocity structure tend to grow into shocks. As 
a result the wind of this type of star is a mixture of rari- 
fied hot material and colder dense regions, with "cool" and 
"hot" gas in dynamical equilibrium. It is common in the 
literature to attribute wind "clumps" to dense formations 
with temperatures much lower than the temperatures of 
the hot gas. These clumps are presumably responsible for 
observed line-profile variability or the formations of DACs 
(discrete absorption components) observed in the optical 
and UV spectra of O stars. The estimated quantity of 
these clumps is of the order 10^ (Howk et al. 2002); their 
size hierarchy varies and appears dependant on spectral 
type and possible structure of local magnetic fields. The 
largest DACs can be seen in the line spectrum of the star 
for ^few days but do not exist longer than a few rotational 
periods of the star (Kapcr ct al. 1997) 

The whole Cyg X-1 system, including the black hole 
secondary, is surrounded by the relatively dense clumpy 
stellar wind. Its large amplitude of variability changes the 
column density of the neutral gas significantly in the di- 
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rection of the secondary. Two types of such variabihty are 
known: orbital, that is connected with the variable col- 
umn density due to orbital movement of the secondary 
in the envelope of the primary component (Brocksopp et 
al. 1999b; Wen et al. 1999); and dips, which usually last 
in the X-ray spectrum for several minutes up to 8 hours 
(Balucihska-Church et al. 2000). Both types of X-ray vari- 
ability are seen during the low/hard state. Dips in the soft 
X-ray spectrum are relatively widely spread. During the 
dip, the spectrum is hardened and an iron K absorption 
edge may be seen, showing that they are due to photo- 
electric absorption (Kitamoto et al. 1984). Outside the 
dips the column density of neutral gas in the direction 
of secondary is ~ 6.21 ± 0.22 x lO^^cm-^ (Schulz et al. 
2002) but can easily reach ~ 1 x lO^'^cm^^ during dips 
(Baluchiiiska-Church et al. 1997). The size of the clumps 
that are responsible for short-lived dips is ~ 10^ cm, i.e. 
relatively small in size. The quantity of clumps on the line 
of sight depends on the geometrical length of the envelope, 
so dips are more frequent at superior conjunction (Wen et 
al. 1999; Baluchihska-Church et al. 2000). 
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